Entanglement enhancement in spatially inhomogeneous many-body systems 
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We investigate the effects of spatial inhomogeneities on the entanglement of modes of strongly 
correlated systems in the framework of small Fermi-Hubbard chains. We find regimes where entan- 
glement is strongly enhanced by the presence of inhomogeneities. This contrasts recent reports of 
entanglement destruction due to inhomogeneities. We further study this phenomenon using concepts 
of Density Functional Theory and, thus, provide a general recipe for the prediction of entanglement 
enhancement in nanostructures. We find enhancement of up to ~ 24%, as compared to impurity-free 
chains. 
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Nanoscale spatial inhomogeneities are ubiquitous in 
strongly correlated systems. Either naturally occurring, 
as in disordered media, or artificially prepared in nanos- 
tructures and ultracold atoms, spatial inhomogeneity can 
influence optical, electrical, magnetic, transport and en- 
tanglement properties [THT7]. Therefore, it is of great 
fundamental and technological interest. 

The impact of spatial inhomogeneities on entangle- 
ment, which is a key ingredient in quantum information 
theory, has been investigated in spin chains, Kondo and 
Hubbard models [T2HT9] . The main goal is the devel- 
opment of solid-state devices for quantum information 
processing [I4j [20] , for which accessing high degrees of 
entanglement, despite the inhomogeneities, is essential. 

For the important case of the inhomogeneous ID 
fermionic Hubbard model, the entanglement of modes 
[2T| [22] or occupation-number entanglement has been 
studied. In this case, the ground-state entanglement be- 
tween two sections (or subsystems) of the chain is quan- 
tified by the von Neumann entropy (see Eq. (2) below). 
The influence of several types of inhomogeneities were 
analyzed: localized impurities and superlattice struc- 
tures [17] . harmonic traps [T7HT9] and disorder [12 j. Al- 
though entanglement depends on the specific inhomoge- 
neous component of the potential, the effect of inhomo- 
geneities was, remarkably, found to be qualitatively the 
same in all cases: entanglement decreases. 

All those studies focused, however, on the entangle- 
ment between a single site — the smallest possible sub- 
system — and the remaining lattice sites. The scaling 
relations from the single-site entanglement to the entan- 
glement of larger subsystems, called block entanglement, 
have been extensively explored from the statistical point 
of view [23-29 . Nevertheless, it is still a fundamental 
open question whether the block entanglement is also de- 
stroyed by spatial inhomogeneities. Furthermore, from 
the view point of realistic applications, it is crucial to 
build a detailed characterization of the regimes, at which 



entanglement is strong regardless of unavoidable inhomo- 
geneities. 

In this paper, we address this issue by investigating the 
block entanglement of small ID Fermi-Hubbard chains, 
where the inhomogeneity is induced by an external po- 
tential, simulating either localized impurities or superlat- 
tice structures. We find that entanglement is not always 
destroyed: there are regimes where it can be even en- 
hanced by the presence of inhomogeneities. We explore 
the physics behind this phenomenon from the perspective 
of Density Functional Theory, and estimate the optimally 
achievable entanglement enhancement under rather gen- 
eral conditions. Applying our approach to several super- 
lattices, we predict entanglement enhancement of up to 
~ 24% as compared to the homogeneous case. 

The ID inhomogeneous Hubbard model is given by [30] 

L L L 

h = -t ^2 4* d jv + U Y1 ™*t% + ( 1 ) 

<ij>,a i i,a 

with on-site repulsion U and repulsive impurities of 
strength V*, which we will both measure in units of the 
tunneling strength t between adjacent sites. We adopt 
balanced populations (N^ = N±), where N = + N± 
is the total number of particles, L the lattice size, and 
n = N/L the filling factor. The ground-state entangle- 
ment between a block of x lattice sites and the remaining 
L — x sites is quantified by the von Neumann entropy 

S x = -Tr [p x log 2 p x ] , (2) 

where p x = Ti{ L _ x y [p] is the reduced density matrix, and 
p, which is calculated by exact Lanczos diagonalization, 
is the ground state's total density matrix. 

To set the stage, let us recall the general properties 
of the entanglement as defined by S x : The latter de- 
pends on the degree of purity of the reduced density ma- 
trix, when it is pure (p x = \ip x )(' l Px\)i S x — 0; an d the 
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FIG. 1. Illustration of the system with 2 impurities (V) and 
the relevant parameters for a block x: impurities (a) sym- 
metrically distributed (2 imp sym); (b) asymmetrically dis- 
tributed (1 imp asym) and (c) at a distance d. 



more mixed p x , the greater S x . That means entangle- 
ment vanishes whenever one of the possible states (in 
the reduced system) has unitary probability, whereas the 
maximum entanglement is achieved when all reduced sys- 
tem's states are equally probable. For the infinite homo- 
geneous Hubbard chain, in the basis of occupation num- 
ber, this condition is fulfilled at half- filling (n = 1) and 
U = [3T] . Thus, for a fixed x and periodic boundary 
conditions, conformal invariance predicts the maximum 
entanglement to be [28j [29] : 



? hom 



(x, L -> oo) = 2 + log 2 [ai 2/3 ]. 



(3) 



Localized impurities. We start by analyzing a system 
with 2 neighboring repulsive impurities of same strength 
V, as Figure 1 illustrates. We use periodic boundary 
conditions, such that the location of the impurities in 
the chain is not important. The relevant parameters are 
the impurities' location within the block x, if x contains 
impurities; or the distance d from the impurities, when 
x does not contain impurities. 

The block entanglement is addressed from three dis- 
tinct perspectives. We consider the entanglement i) be- 
tween impurities and the remaining sites, the so-called 
impurity entanglement [I5j [16], and ii) between bipar- 
titions containing both impurity and non-impurity sites. 
We also study Hi) the average over all possible bipar- 
titions, which quantifies the overall impact of inhomo- 
geneities (incorporating i and ii), particularly important 
for the treatment of disorder. 

Figure 2 shows the results for the single-site entangle- 
ment as a function of V. We find that only the impu- 
rity entanglement (case i) vanishes for V » t, as the 
reduced density matrix of the impurity sites is pure in 
this limit (p x = pimp — |0)(0|). For the average over all 
possible bipartitions (case Hi), we find that the entangle- 
ment decreases with V, but saturates at a finite value, 
what is consistent with previous observations [I2j [T7HT9] . 
Our results in Fig. 2 reveal that the main contribution 
to this decrease comes from the impurity entanglement. 
Surprising, though, is the behavior in the situation when 
blocks do not coincide with impurity sites (case it), as 
can be seen in more detail in the inset of Fig. 2. In this 
case, the degree of entanglement is not reduced. On the 
contrary, it is even enhanced by V > 0. We study this in- 
teresting entanglement enhancement based on a Density 
Functional Theory (DFT) [32H34] point of view. 
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FIG. 2. (Color online) Single-site entanglement (x — 1) as 
a function of the impurity strength V/t for the impurity site 
(imp), for single sites at distance d from the impurities, and 
the average over all possible bipartitions. Here L — 10, N — 
6, U = 4t, and we adopt periodic boundary conditions. Inset: 
enlarged view of entanglement enhancement as compared to 
the homogeneous case (dashed line). 



Within DFT, all observables are functionals of the 
single-particle density n(r) of the many-body interact- 
ing system. Unfortunately, the explicit form of the den- 
sity functional of a given physical entity is, in general, 
unknown. However, for homogeneous Hubbard chains, 
where all the sites are equivalent (n(r) = n), the den- 
sity functional for the single-site entanglement is exactly 
obtained from a Bethe-Ansatz solution [31]. One finds 
a broad regime of densities (n < 0.8 for U > 0) where 
the entanglement monotonously grows with density. For 
inhomogeneous systems, which have by definition site- 
dependent densities (n(r) = {rii}), it has been proven 
that the homogeneous functional is a good approxima- 
tion for the single-site entanglement, within a local den- 
sity approximation (LDA) replacing n — >• rii [I2j [T71 [35] . 
Consequently, the inhomogeneous single-site entangle- 
ment varies with density as in the homogeneous case [31] : 
monotonically increasing with density for rii < 0.8. 

This is precisely the mechanism behind the entangle- 
ment enhancement observed in Fig. 2: the particles are 
repelled from the impurities, inducing a larger effective 
density at the remaining L — 2 sites (n eff = N / (L — 2) > 
n), which thus leads to higher entanglement, by virtue 
of the above DFT argument [36 . This can be verified 
at strong impurity strengths, where the impurities repre- 
sent physical boundaries and therefore mimic an effective 
system with open boundary conditions. In Figure 3a, we 
compare the density profiles of such effective chains to 
our original model with impurity strength V = St: the 
very good agreement confirms our interpretation. 

There are some additional features in Fig. 2 which ne- 
cessitate some refinement of our above LDA argument. 
To start with, the entanglement of blocks d = and d = 1 
differs by - 6% at V = St (inset of Fig. 2), while the 
difference between their densities is negligible (Fig. 3a). 
Furthermore, block d = 2 has the smallest density among 
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FIG. 3. Density profile of the system with 2 impurities (see 
Fig. 1), at V = St: (a) The blocks considered in Fig. 2 are 
identified by d, and the local densities (circles) are compared 
to the densities in an effective chain of reduced length L = 8, 
V — (squares), with otherwise identical parameters as in 
Fig. 2, and subject to open boundary conditions, (b) Pre- 
dictions of the optimal entanglement enhancement for several 
sizes x as highlighted by the coloured shapes, according to 
our interface-density amended LDA approach. 



the non- impurity blocks (Fig. 3a), but it does not exhibit 
the weakest entanglement (Fig. 2). The most curious 
feature occurs for the block at d = 0: its entanglement 
behaves as the other d > blocks only for very small 
values of V and then decreases considerably (see inset of 
Fig. 2). All these features of the ^-dependence of entan- 
glement ought to be encoded, at least in principle, in the 
exact density functional for the entanglement, which is a 
highly nontrivial function of {n^}. 

The simplest possible functional beyond LDA is to as- 
sume that the block entanglement is not only determined 
by the block density itself, but is also affected by the den- 
sity at its border sites, i.e., at the outmost sites of the 
complementing block L — x, in accord with the area law 
concept [37]. In our case, the boundary surface is fixed, 
but weighted by the local particle density. Thus, incorpo- 
rating the boundary density into our above relation be- 
tween entanglement and local density, we conclude that 
the larger the interface density (up to ~ 0.8), the larger 
its contribution to block entanglement. 

This simple, interface-density amended LDA approach 
is already enough to recover qualitatively all features of 
Fig. 2. In particular, it explains the special behavior 
at d = 0: one of the border sites contribution is sup- 
pressed (77,4 ~ 0, see Fig. 3a), such that the entangle- 
ment saturates at a smaller value compared to all the 
other d > cases, in which both border sites (with non- 
vanishing densities) contribute to entanglement. 

Extending this analysis to x > 1, we can now in- 
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FIG. 4. (Color online) Block entanglement as a function of the 
impurity strength V/t for (a) x = 2, (b) x = 3 and (c) x = 4, 
at different distances d and symmetrically (sym) or asymmet- 
rically (asym) distributed impurities (also see Fig. 1). (d) 
Average over all possible bipartitions. L = 10, N = 6, U — At 
and periodic boundary conditions. 



fer which block exhibits the maximum entanglement en- 
hancement for a given block size x. Figure 3b indicates 
the blocks with optimal entanglement enhancement for 
our case with 2 impurities. Figure 4 confirms this expec- 
tation: The blocks with optimal entanglement enhance- 
ment are those with the highest interface density. Quali- 
tatively all the features observed at x = 1 are seen again 
for x > 1 (except that, for x > 2, as the block is larger 
than the number of impurities available, the impurity- 
entanglement scenario is not defined), but the precise 
scaling relation from x = 1 to larger x remains non-trivial 
and unknown. The behavior of the average block entan- 
glement (perspective Hi), however, is similar for all x, 
as Figure 4d shows. This quantity contains less details, 
consequently, simpler scaling properties. 

Superlattice structures. We now apply the general fea- 
tures of optimal entanglement enhancement to superlat- 
tices (SL) pQ. Here, we consider SL structures as defined 
by a periodic modulation of V*. From our previous anal- 
ysis, one expects substantial entanglement enhancement 
in superlattices, whenever block and remainder are com- 
posed by border sites with V = and inner sites with 
V > (similar to the x = 4 case, Fig. 3b). 

Figure 5 shows the SL entanglement for blocks with 
this specific configuration in several distinct SL struc- 
tures. We find entanglement enhancement for all cases 
in which the effective density (n eff = N/(L — I), I 
the total number of impurities) is larger than n (up to 
n eff ^ 0.8), consistently with the localized impurities re- 
sults. Remarkable though is the much higher relative 
enhancement: entanglement up to 24% larger than in 
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FIG. 5. Entanglement for several superlattice structures for 
x — 3 and x — 4: SL^ with period [a, 6, c, c2], where up- 
per indices are sequences of sites (layers) with V% = V, while 
lower indices are layers with V% = (for instance, the unit cell 
for SL^ and SL^ are [V0VV00] and [V00], respectively). In 
all cases the lattice parameters are U = 4£, V = St and pe- 
riodic boundary conditions. The entanglement enhancement 
is relative to the homogeneous chain (V = for all sites) of 
same length L: The upper data set (*SL^>) is obtained from 
DMRG calculations for L = 36, while all others from Lanczos 
diagonalization for L — 12. 



the impurity- free system is observed. 

The absolute entanglement values in SL are, in some 
cases, even higher than the maximum homogeneous en- 
tanglement one can achieve by properly adjusting the 
parameters, as predicted by Eq. (3) (Sy^(x = 3) = 3.06 
and S^£(x = 4) = 3.33). This certainly encourages 
the manipulation of nanostructures towards progress in 
quantum information devices. 

We have also considered a few cases of SL entangle- 
ment for larger periodic chains (up to L = 36), by using 
Density Matrix Renormalization Group (DMRG) tech- 
niques [38]. We find that the impact of L on the entan- 
glement enhancement is weak, as can be seen in Fig. 5, 
by comparing SL3> (for L = 12) and *SL3> (for L = 36). 

In summary, our results demonstrate that spatial in- 
homogeneities can actually act in favor of entanglement 
and, accordingly, could be cleverly engineered in solid- 
state devices for optimal quantum information processes. 

We thank Dominik Horndlein for providing us with his 
DMRG code. VF is indebted to CAPES (4101-09-0) and 
to the Forderung evaluierter Forschungsprojekte of the 
University Freiburg. 
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